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Abstract 



The continuous block spin (Wilson) renormalization group equation gov- 
erning the scale dependence of the action is constructed for theories con- 
taining scalars and fermions. A locally approximated form of this equation 
detailing the structure of a generalized effective potential is numerically an- 
alyzed. The role of the irrelevant operators in the nonperturbative renor- 
malization group running is elucidated and a comparison with the 1-loop 
perturbative results is drawn. Focusing on the spontaneously broken phase 
of a model possessing a discrete symmetry forbidding an explicit fermion 
mass term, mass bounds on both the scalar and fermion degrees of freedom 
are established. The effect of the generalized Yukawa coupling on the scalar 
mass upper bound is emphasized. 
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1 Introduction 



A basic, generally unstated, tenet used when describing elementary particle 
interactions is that the physics at one distance scale follows uniquely from 
the dynamics specified on a more finely grained distance scale. As such, it 
is crucial to understand the short distance dynamics of the standard model 
of elementary particle interactions. Since the non-Abelian gauge interac- 
tions are asymptotically free, their short distance behavior is amenable to a 
renormalization group improved perturbative analysis. On the other hand, 
the scalar self interactions and the scalar-fermion Yukawa interactions are 
characterized by couplings which are not asymptotically free. Consequently 
any conclusions drawn from a perturbative analysis of their short distance 
structure is at best indicative and require confirmation using nonperturbative 
techniques. 

For instance, a study of the one-loop effective potential in the self coupled 
A0 4 scalar model produces a Landau singularity' 1 ' which has in turn been 
used to place a perturbative triviality bound' 2 ' on the scalar mass. The 
qualitative nature of these conclusions has been largely substantiated by 
various nonperturbative studies. ' 3 ~4 Here it is found that a nontrivial self 
coupled scalar model can only be consistently defined as a cutoff theory; that 
is, as an effective description of the dynamics below some momentum scale 
A. The only consistent A — > oo limit of this model is one characterized by 
a vanishing renormalized coupling. This free field behavior is a reflection of 
the Landau singularity already observed in renormalization group improved 
perturbation theory. 

Within the perturbative framework, given the value of the four scalar 
self coupling at the scale A, A (A), which in general is such that A (A) > 1, 
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then the coupling at smaller momentum scales follows from the one-loop 
renormalization group equation as 

1 1 3 A 

A(m^) ~~ A(A) + 11 ~nv H ' 

where m# is the scalar mass. In general, this equation is inadequate as it 
is being applied in a range outside the domain of validity of perturbation 
theory. If however, ran « A, then the low energy value A(m#) is insen- 
sitive to the initial value A(A) due to the renormalization group trajectory 
focusing effects of the trivial (Gaussian) infrared fixed point for A. On the 
other hand, when discussing upper bounds on the scalar mass, one must con- 
sider values of itlh which are not much smaller than A itself. For instance 
for A ~ 1 Tev, the upper bound on m# using perturbation theory' 2 ' is found 
to be m H ~ 800 GeV. Certainly these quantitative results require substan- 
tiation in a nonperturbative framework. In fact it turns out that various 
such estimates' 3-4 ! have been obtained yielding a maximum scalar mass of 
m H ~ 650 — 800 GeV again for A ~ 1 TeV. This rough quantitative agree- 
ment of these nonperturbative predictions with the perturbative estimate is 
certainly not a priori anticipated. In addition to the fact that the scalar self 
coupling is sufficiently large to render the 1-loop perturbative result suspect, 
it also follows that such a coupling induces additional irrelevant operators 
whose effects are not necessarily insignificant for scales not too far removed 
from A. In general, the coefficients of these operators are suppressed by pow- 
ers of m H /A. While for m H << A such terms are clearly quite small, this 
need not be the case as m H approaches A. 

Of the various nonperturbative studies which have been performed for 
these scalar self coupled systems, most have employed lattice Monte Carlo 
simulations. ^ A somewhat different approach' 4 ] to the study of these systems 
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employs the continuous Wilson renormalization group. ' 5 ~ 7 ' In this case, one 
specifies the model at the ultraviolet (UV) cutoff scale A and, by demanding 
the cutoff independence of the renormalized Green functions, one constructs 
a functional differential equation for the action at the longer distance scale. 
This action contains not only the relevant and marginal operators but also 
the irrelevant operators with the coefficient of each operator fixed in terms 
of the initial action defined at scale A. Moreover as the construction is 
independent of the strength of the initial coupling, it constitutes a nonper- 
turbative procedure. Using a local action approximation' 8 ' to these exact 
functional equations, Hasenfratz and NagerW extracted an upper bound for 
the scalar mass consistent with the other nonperturbative estimates. Re- 
cently, the sensitivity of this bound to the choice of action at scale A has 
also been investigated.^ The resulting maximum scalar mass was found to 
be insensitive to the initial action. 

The purpose of this paper is to extend the Wilson renormalization group 
approach to the construction of renormalized effective Lagrangians in theories 
containing fermions in addition to scalars. Once again, a consistent model 
of this sort requires the introduction of an ultraviolet cutoff A. Provided 
that a huge hierarchy exists between A and the scalar and fermion masses 
(generally accomplished by a fine tuning of parameters), then the infrared 
stable quasi fixed point' 10 ' of the Yukawa and scalar self couplings again fo- 
cuses the renormalization group trajectories to a low energy value basically 
insensitive to the exact nature of the short distance couplings. Moreover the 
position of the infrared quasi fixed point can be determined using the per- 
turbative 1-loop renormalization group equations for these couplings. This is 
precisely the situation which has been exploited in the top quark condensate 
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models, t 11 ' 

On the other hand, for scalar and fermion masses not too far removed 
from the UV cutoff, some nonperturbative tools must be employed. In the 
next section, we use the continuous Wilson renormalization group to con- 
struct a functional differential equation for the action functional which can 
serve as just such a nonperturbative procedure. To obtain this equation, one 
varies the UV cutoff and demands the resultant action functional reproduces 
the same physics (i.e. renormalized Green functions) on all momentum scales 
less than the new (lower) cutoff. That is, the change in cutoff is compensated 
by changing the coupling constants of a complete set of local operators. This 
set includes irrelevant operators in addition to the relevant and marginal 
operators. Thus the various action functionals satisfying this equation cor- 
respond to points on a particular (Wilson) renormalization group flow each 
producing the same physics on momentum scales less than the UV cutoff 
used in its definition. 

This exact Wilson renormalization group equation is tantamount to an 
infinite number of coupled consistency equations for the various coupling 
constants of the complete set of operators. As such, its analysis and extrac- 
tion of physical results from it requires additional approximations. In section 
3, we employ a local action approximation in terms of which the functional 
differential equation reduces to a partial differential equation for a general- 
ized effective potential in a model containing one scalar and one fermion. 
In the limit of small scalar self coupling and Yukawa coupling at the UV 
cutoff scale, this equation reduces to the familiar 1-loop perturbative renor- 
malization group equations for these couplings. For larger couplings at the 
UV cutoff scale, however, the equation for the generalized effective poten- 
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tial requires a numerical evaluation. We perform such an analysis under the 
further assumption that only terms up to bilinear in the fermion fields are 
retained in the generalized potential. The resulting numerical solutions are 
then analyzed for various choices of the initial Yukawa and scalar self cou- 
plings and for different initial cutoffs. In particular, the role of the induced 
irrelevant operators in driving the theory toward the trivial (Gaussian) fixed 
point is emphasized. The exact numerical solutions are also contrasted with 
the analytic solutions of the one loop perturbative equations extrapolated 
well outside their domain of validity. 

By examining the numerical solutions resulting from different values of 
the scalar self coupling and Yukawa coupling chosen at various cutoff scales 
and resulting in nontrivial scalar vacuum expectation values, we obtain in 
section 4 the trivialty and vacuum stability bounds on both the scalar and 
fermion masses. Our results indicate that the scalar mass bound increases in 
the presence of the Yukawa coupling. The amount of this increase depends 
on the strength of the initial scalar self coupling; larger self couplings yielding 
smaller percentage increases. For the range of couplings we considered, the 
percent increase in the scalar mass bound varied from roughly 15-20% to 6- 
10%. Similar conclusions have also recently been reached from lattice Monte 
Carlo simulations of scalar- Yukawa theories. ^ In addition, using the values 
of the mass upper bounds as a function of the cutoff, we determine the allowed 
ranges of scalar and fermion masses for a specific choice of the cutoff. 
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2 Wilson Renormalization Group Equation 

Consider a general Euclidean quantum field theory containing the fields fij. 
The subscript % delineates bosonic or fermionic fields, as well as labels any 
Lorentz and internal symmetry indices, e.g. Jl, G {ip,^,^, A^, ...}. The 
theory is regulated so that the contributions to the Schwinger functions (the 
Euclidean space Green functions) of the degrees of freedom with momentum 
above A are strongly damped, and it is renormalized so that the Schwinger 
functions are independent of the cutoff A. Thus the generating functional for 
the renormalized Schwinger functions, Z[J], is cutoff independent: j^Z[J\ = 
0. In general, Z[J] may be represented by the path integral 



Here S[Q; A] is the Euclidean action, assumed known at scale A where the 
coupling constants are specified while the sources Jj are of complimentary 
character to the fields Qi and have support only on momentum scales below 
A. More specifically, the fields depend on the chosen regulation scheme, 
Qi = Qi(p; A), while the action has the form 



with given coupling constants g(p;A). The detailed question of regulator 
dependence of S will be addressed in future work. For simplicity, in the 
following, we consistently employ a hard cutoff regulator in momentum space. 
Thus the degrees of freedom with momentum above A are taken to be zero 
and we can write Q(p; A) = 6>(A — |p|)Q(p). The action then becomes 




(2.1) 




9h-i n (Pi, ■ ■ ■ ,Pn, A) ^n(-Pi; A) • ■■n in (-p n ] A) , (2.2) 
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SP;A] = E^/^j--^P^ 4 fe + --+p»)«(A-M) 

• • • 0(A - \PnD9h-in (pi, • • • , A)^u (-Pi) • • • O in (-p n ) , 

(2.3) 

while the normalized functional measure in Eq. ( |2.1| ) contains only those 
degrees of freedom with momentum less than A: [dQ] = I\\ P \<Ad^{p)- 

To derive the continuous Wilson renormalization group equation, we al- 
low the cutoff to vary and consider the response of the theory subject to the 
constraint that the Schwinger functions remain unaltered. Parametrizing the 
cutoff as A(i) = e~'A, the dependence of the action on A(t), or equivalently 
on t, can be determined by integrating out the degrees of freedom with mo- 
mentum lying between the differentially close cutoff values, A(t = 0) = A and 
A(5t) = e~ St A = A — 5tA. Towards this end, we expand the action 5* around 
the degrees of freedom with momentum below A(5t) retaining terms up to 
linear in 5t. This requires retaining only terms no higher than quadratic 
in the fields with momentum in the differential shell e~ St A < \p\ < A since 
each momentum integral over the momentum shell introduces a differential 
factor of 5t (save for the one momentum integration over the overall energy- 
momentum conserving delta function). This Gaussian integral over the fields 
with momentum in the shell can be performed yielding an action correspond- 
ing to the cutoff A(St), which can then be used in the path integral for the 
remaining lower momentum fields. So doing we find 

z[j] = f idn]e- s[Q - Am+ f $? J ^- P ) {2A) 

J\p\<A(8t) 

The new action defined at the lower cutoff A(5t), S[Q; A(6t)], is thus obtained 

as 
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If d 4 p (M(p,-pY 



S[Q;A(5t)] = 5[n;A(0)]| + -/ T^strln 

I J shell \^Tl) 

r d 4 p d 4 q F . 5S[Q 

'shell (27r) 4 (27r) 4 ^ 1 ~SQ~(p) 



_1 f d^tfq_ Fi 5S[Q; A] ft?[n;A] 
2 y ahe jj (2nY (2nY 1 ; (Jafp) lj {P ' q) 6£lj(q) 

(2.5) 



where Fi is the Grassmann number (even for bosons, odd for fermions) of 
the field Qi and the vertical lines appearing on the right hand side indicate 
that each such term is to be evaluated setting the fields having momentum 
greater than e~ 5t A to zero. We have further defined the inverse propagator, 
Mij(p, q), as the second order expansion coefficient of the action so that 



M-(p <,) = (-)<* ^ [Q;A] 



(2.6) 



where Ay is unity if both the field is an even element of the Grassmann 
algebra and the field Qj is an odd element of the Grassmann algebra and 
vanishes otherwise. More explicitly, M has the form 

(M BB M BF \ 
M — I , 2.7 

\M FB Mpp J 

where M BB (Mff) is the submatrix of bosonic character obtained by differen- 
tiating S with respect to two bosonic (fermionic) fields, while the submatrices 
M BF and M FB are of fermionic character and result from one bosonic and 
one fermionic field derivative. The superdeterminant of M is defined in terms 
of the supertrace of the above component matrices (see the Appendix for a 
review of manipulations with Grassmann valued matrices) yielding 

, , , detM BB 
sdetM = - j , 2.8 
detAW ' v ; 
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with 

N FF = M FF - M fb MbIM bf , (2.9) 

and hence the result 

str InM = tr \nM BB - tr lnN FF . (2.10) 

The action S[f2; A] | is obtained from S^fljA] by setting all modes with mo- 
mentum above A(5t) to zero and thus can be expanded as 

^;A]| = ? -L/^L...^( 2 .)V(p 1 + ... + p„) 

9{A(5t) - \ Pl \) ■ ■ ■ 6(A(6t) - \ Pn \)g iv .. in (pi, ■ ■ ■ ,Pn, A) 

n il (- Pl )...n in (- Pn ). (2.11) 

(Note that the coupling constants are still those originally defined at A.) The 
new action appearing on the left hand side of Eq. (|2.5| ), on the other hand, 
involves the lower momentum fields as well as the coupling constants defined 
at the lower cutoff and can be written as A(5t) 

mm) = e^/^-^(2^v + - + p.) 

9(A(6t) - M) • • -9(A(6t) - \ Pn \)g n .., n (Pu A(<5t)) 
Q H (- Pl )...Q ln (-p n ). (2.12) 

Their difference just results in the variation of the coupling constants when 
the cutoff is lowered so that 

S[Q;A(6t)} - ^A]|=E^/^---^(2vr) 4 5V + ---+ P „) 



6(e- dt A - | Pl |) • • • 9(e~ dt A - \p n \) [g( Pl , ...,p n ; e~ dt A) 

-g(pi,..., P n \ A)] n fl (~Pi) ■ ■ ■ fli„ (-p n ) 
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6(A-\ Pl \)...6(A-\ Pn \) 
d 1 



(2.13) 



where the overall factor of 5£ allows t to be set to zero inside the integrand in 
the last line on the right hand side (we consistently retain only terms through 
5t). 

In order to obtain a differential equation for the action functional S it- 
self, the A— derivative must be brought past the regulated fields (i.e. the 
^-functions) . This can be achieved by rescaling all the momentum vectors to 
unit vectors 

p^ = Ap^ (2.14) 

so that 9 (A — \p\) = 0(1 — \p\) and by rescaling all the fields to dimensionless 
ones as 

^i(p) = C£(A)aG3) , (2.15) 

(i.e. for the regulated fields, f^(p;A) = Cq 1 (A)&i(p; 1)). It follows that the 
coupling constant variation now includes a momentum dependence since 

d d n d 

A—g( P] A) = A—g(Ap; A) - ^—^Ap; A) . (2.16) 

Substituting these variable changes into Eq. ( |2.13| ), and recalling that A(t) = 
e~ l A so that a differential logarithmic change in A corresponds to a negative 
differential change in t, we obtain 
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A— S[n;A(5t)] 



d 
dt 



S[Cl;t] 



6t 



i{5[fi;A(5t)]-5[fi;A]|} 
d A p 



4-E 

n, 



(27T) 



&\ 5 



S[Cl;A] 



(2.17) 



where S is the same action as S, but expressed in terms of the rescaled 
quantities as 



S[£l;t] = S[Q;A(t)] 



^n!7 (2vr^ 



(i 4 Pi d 4 j3 n 



(2vr) 4 5 4 (p 1 + ---+K, 



(2vr) 4 

0(l-|p 1 |)...0(l-|p n |) 

<7ii-i„ (pi, • • • , Pn] A)Q h (-pi) ■■■Q in (-Pn) , (2.18) 



with the dimensionless coupling constants defined according to 



g il ... in (p 1 ,...,p n ;A) = 



A 4 "- 4 I1^.(A) 



gi x ...i n {Ap-i, ...,Ap n ;A) 



(2.19) 



In addition, the anomalous dimension, 7^, of each field is defined by 

d 



A—\nCa i (A) = 4-dn i +7n i 



(2.20) 



with dn t the naive field dimension; e.g. dfi. = 1 for bosons, and da f =3/2 
for fermions. The prime on the momentum derivative in Eq. ( |2.17| ) indicates 
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that it does not act on the overall energy-momentum conservation Dirac 
delta function or ^-functions but only on the dimensionless couplings in the 
expression for S. 

Finally, substituting Eq. ( |2.5|) with M expressed in terms of the rescaled 
momenta and fields into Eq. ( |2.17 ), we obtain the continuous Wilson renor- 
malization group equation for the action: 



d_ 

It 



S[£l;t] 



str lnM(p, — p) 



11 r d 4 p 

2 6t Je-s*<\ P \<i (2vr) 4 " 
1 1 f d 4 p d A q Fi 5S[fl;A] 

-"<HM<i (2 7 r) 4 (27r) 4(_ ' ) ' 

SS[n ; A] 



2 St 



+ 



4-y / **- 

tr4l<l (2vr) 4 



dn t - In, + p 



S 



dp^ J Stti(p) 



A] . 



(2.21) 



Here and in what follows we drop the carat notation on all quantities and 
work only with dimensionless fields and unit momentum vectors. 
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3 The Local Action Approximation 



The Wilson renormalization group equation ( 2.21| ) is an exact equation de- 



termining the action at scale A(t) given its initial short- distance form at 
A = A(0). Unfortunately its complexity is such that it can only be solved 
approximately at the present time. Towards this end, we ignore all inter- 
actions involving derivative couplings and furthermore set the wavefunction 
renormalizations to unity. This is clearly a drastic approximation but does 
not leave the equation completely devoid of content. Actually, in the small 
coupling constant regime, this approximation reproduces the one-loop per- 
turbative renormalization group equations with anomalous dimensions ne- 
glected. The locally approximated action, S[fl; A], takes the general form 



where the coefficients v\J.. in are momentum independent and it is understood 
that the Euclidean momentum integrals are over unit vectors. 

When expanding the fields about space-time independent values, fli(x) = 
Qi, (or in momentum space Qi(p) = (2ir) 4l 5' i (p)Qi), it is convenient to intro- 
duce the "generalized potential" U (fl; A) as 




1 r d 4 pi d 4 p 2 
2J (2tt) 4 (27r) z 

+ E; 



1 (2n) 4 5\p 1 +p 2 )n il (-p 1 )K tll2 (p 1 )n t2 (-p 2 ) 




n 



toii(-Pl)~-£Un{-Pn), (3.i) 



17(0; A) 



n 



—v {n) n ■■■q 



(3.2) 
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so that, when evaluated at constant fields, and using Kij(0) = 0, the action 
simply takes the form 

S[Q = constant; A] = (2vr) A 5 A (0)U(Q; A) . (3.3) 

Thus for constant fields, the Wilson renormalization group equation ( |2.21| ) 
in this approximation reduces to a differential equation for this generalized 
potential given by 



9 rr/n x 1 1 f d P 

—U(n,t) = -— / — -V-strln 

dt v ' ; 2 8t Ji-st<\p\<i 2vr 4 



+ 



4 - d nMi 



d 



U(Q,t). (3.4) 



In obtaining this result, we have used the fact that the "tree" graph contri- 
butions to Eq. ( |2.21| ) (the penultimate term on the right hand side) vanish in 
the constant field expansion since it sets the "propagator", M -1 , momentum 
to zero, which is far from the momentum shell 1 — 8t < \p\ < 1. In addi- 
tion, the approximation of neglecting the anomalous dimensions, 7^. = 0, is 
consistent with the assumption of no wavefunction renormalization. 

In order to proceed further, a particular choice of fields and action must 
be made. For simplicity, we will study the field theory of one real scalar 
field if and one Dirac spinor fermion field ip and its conjugate x/j. The local 
approximation action then has the form 

S[ip, ip, if), t] = — - J d^xipd 2 ^ + J d 4 xi])-y ■ dip 



+Eii^ w w/^ , (#r. ( 3 - 5 ) 



l,m 

where for the possible fermion interaction terms we have assumed only prod- 
ucts involving scalar bilinears enter. Thus only functions of iptp occur in U. 
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If the discrete 75-symmetry, 



i>' = -^75 



(3.6) 



is also imposed, the interaction terms involve monomials with I + m equal an 
even interger and so explicit fermion mass terms are forbidden. Hence the 
generalized potential is simply 



Lm 



(3.7) 



and the kinetic energy matrix in {(p,ip,ip} space is 

^ 2 



(P 



K{p) 



■ T 

—vy ■ p 



(3.8) 




V — ry ■ p 

Here the Hermitian Euclidean space 7-matrices are defined to satisfy 
= 26^, with 75 = 7i727374- 
Inserting these expressions in the renormalization group equation ( |3.1| ) 
and performing some straightforward, but lengthy, algebra we find the locally 
approximated form of the continuous Wilson renormalization group equation 
for the generalized potential given by 





Of 



U(ip,a;t) = AU(ip,a;t)-ipU (p (ip,a;t)-3aU a (ip,a;t) 

+ — In (1 + U^b, v, t)) - In (l + U*{tp, a; t 



16tt 

+ In (1 + £(</?, a; t)) 



(3.9) 



16 



where we have introduced the a field as a = ipip and have defined the com- 
bination 

£ = 2^fer £/™ - E/™ I , 1 „ 1 . (3.10) 



^ Ua<p {i + u v J u * a 



1 + ^ 2 

In addition, we have denoted the derivatives of the potential with respect 
to the field variables by subscripts on the generalized potential, so that, for 
example, = -^q^U. Finally, we shall assume that the model at the UV 
cutoff scale is simply given by the canonical kinetic matrix K of Eq. ( |3.8| ) 
and the initial generalized potential function 

U(cp, a; t = 0) = ^/i 2 (0V 2 + ^\(0)<p* + g(0)<pa , (3.11) 

which is invariant under the discrete 75-symmetry of Eq. ( |3.6|) . Note again 
that this symmetry forbids the appearance of an explicit fermion mass term 
or equivalently a term linear in a. 

We seek solutions for the generalized potential in which this discrete 75- 
symmetry is spontaneously broken. This is achieved by choosing the pa- 
rameters of the model at the cutoff scale A such that the vacuum energy 
is minimized by the scalar field if acquiring a nonzero vacuum expectation 
value. At the same time, we assume that no fermion bound states or con- 
densates form. As such, the vacuum energy is identified with the potential 

V(ip; t) = U{ip, a = 0; t) . (3.12) 

That is, the vacuum is characterized by (p 7^ but with a = 0. 

In order to proceed with the numerical solution of the partial differential 
equation, we make the further truncation of retaining terms in the generalized 
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potential only through those linear in a (i.e. ipip) so that 

U(ip,a',t) = V{ip;t)+aG(ip',t). (3.13) 
Here V(ip;t) is the effective potential function defined in Eq. Q3.12| ), while 



G((p; t) is a generalized Yukawa coupling. This approximation, made for 
technical reasons only, neglects those irrelevant operators containing four or 
more fermion fields. Introducing the function 

F(wt)=V v (wt) (3.14) 

and recalling that a vanishes in the vacuum, then the partial differential 
equation for U reduces to the two coupled partial differential equations for 
F and G given by 

dF _ op Tp , 1 F W 1 GG V 

dt * v 16vr 2 1 + F v 2vr 2 1 + G 2 

dt - G tt *** + \fti*\+F v 8^(1 + G 2 )(1 + F,)- {6Ab) 

From the functions F and G, the running scalar self coupling, X(t), and 
running Yukawa coupling, g(t), are defined as 

A(t) = 

g(t) = G v (0;t) (3.16) 

and are plotted in Figures 1-2 for various choices of their initial values at the 
UV cutoff A (t = 0). 



1 Since ipip = a is a nontrivial even element of the Grassmann algebra and since no 
fermion condensate forms, then the general expansion for U terminates after terms of 
order (ipip) 41 = <r 4 . 
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Figure 1: Scalar self coupling X(t) as a function of t for various initial Yukawa 
couplings and cutoff to scalar vacuum expectation values. 



Figure 2: Yukawa coupling g{t) as a function of t for various initial scalar 
self couplings and cutoff to scalar vacuum expectation values. 
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Using Eqs. ( |3.11| ), ( |3.15| ) and ( |3.16| ), it is readily seen that g{t) decreases 
as t increases. On the other hand, A(t) can either initially increase or decrease 
depending on whether g(0) is greater than or less than the value, g c (0), which 
is a function of A(0) and fi 2 (0) and produces a zero initial A slope: ^A(O) = 0. 
This value is readily computed as 

g 2 M = ~ — ~^~57 — \ • (3-17) 

ycV ; 21 + /i 2 (0) V ; 

For g(0) > (<)g c (0), \(t) initially increases (decreases) as is readily observed 
in Fig. 1. In either case, the effect of the initially large (nonperturbative) A(0) 
is such as to induce non-negligible higher dimensional (irrelevant) operators 
in the action and consequently to drive the generalized effective potential 
away from the surface spanned by the relevant and marginal operators only 
as was initially the case (cf. Eq. (|3.11|) ). These induced operators play an 
important role in the dynamics for small to moderate t values. 

As a measure of their importance, we display in Fig. 3 the evolution of 
the dimension 6 operator coupling 

x(t) = ^e- 2t F wwip (0;t). (3.18) 

Here we have included the explicit scale factor of e~ 2t ~ accompanying 
the dimension 6 operator in the definition of the coupling rather than in 
the operator itself. As seen in Fig. 3, this induced coupling rises rapidly 
to a peak value before beginning to tail off as t is increased. For small to 
moderate t values, this coupling can be quite large and is instrumental in 
helping to drive the system towards the Gaussian fixed point at the origin 
of coupling constant space. That is, the effect of this and the other higher 
dimensional operators is to dampen the \{t) coupling so that eventually the 
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Figure 3: The t-evolution of the coupling \(t) of a dimension 6 operator 
for various ratios of the cutoff to scalar vacuum expectation value when the 
initial couplings are A(0) = 5 and g(0) = 1. 

theory evolves into a region where it can accurately be described using a 
perburbative expansion for the generalized effective potential. 

It is clear from Fig. 1 that there is a range of t values greater than a 
value i t* for which A(t) can be well approximated by a linear fuction of t. 
Over this range, one can therefore write 

U(<p, a; t) = U(ip, a; t*) + (t - f)^U(<p, a; t*) . (3.19) 

Moreover, the coefficient of (t — t*) can be directly gleaned using the locally 
approximated Wilson renormalization group equation ( |3.9D for the general- 
ized effective potential. So doing, one readily computes for this range of t 

2 A more precise definition of t* will be provided in section 4. 
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values that 



A(t) = \(t*) + (t-t*y 



Lpipipipip 



(0;**; 



108A 2 (f 



+ (t-f 



6tt 2 



96tt 2 [l + i^(0;t*) (1 + F v (0; t*)) s 



3g\t*)-2g(f)G (pipip (0;t*) 



(3.20) 



This same result for A(i) can also be secured by computing the generalized 
effective potential for this range of t > t* values using the 1-loop perturbative 
result and then retaining terms linear in t — t*. The 1-loop contribution to 
the generalized effective potential is given by 



U((p,a;t) = U(cp,a;f) 

l r d 4 P 



+- 



-C*-**)<|p|<i (2tt)' 



str In 



d 2 



13 ' dQjdQi 



(3.21) 



Once again if we retain only the terms through linear in a and decompose 
U((p,a;t) as in Eq. ( J3.13Q , we isolate 



-2 



d A p 

-(*-**)<|p|<i (27r) 4 
d 4 p 



In 



p 2 + V^{ip;t* 



-(*-**)<|p|<i (2-k) 



rln 



Gir-.i) = G^-.n+^G^t*) 



p' + G\ v -e) 



d A p 



(3.22) 



-G(wt*)Gl(tp;t*) 



d 4 p 



e-(*-**)<b|<i (2tt)V + V w (ip;t*) 
1 



■(*-**)<bl<i (2tt) 4 {p 2 + G 2 (v;t*))(p 2 + V^ip-t*))' 



(3.23) 

The momentum integrations are standard and can be staightforwardly per- 
formed. Retaining only terms linear in (t—t*), which constitutes a reasonable 
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approximation for the t range under consideration (0 << t — t* < 1), the 



expression for A(t) is found to be identical with that of Eq. ( 3.20 ). It is im- 
portant to note that the 1-loop perturbative result includes the effects of the 
higher dimensional operators. Their presence is crucial for allowing a smooth 
joining with the numerically generated solution of the locally approximated 
Wilson renormalization group equation. 

Of course, as t continues to increase to even larger values, the effects 
of all the higher dimensional operators will eventually decouple and can be 
completely ignored. (They are truly irrelevant!) In addition, the marginal 
operators will have couplings well in the perturbative region. The onset of 
this behavior is more rapidly achieved the larger the value of the initial cutoff 
A. For this range of t values, the evolution of the model can be accurately 
parametrized by summing the leading logorithms of the 1- loop graphs thus 
obtaining the 1-loop perturbative renormalization group running for X(t). 

To perform the perturbative analysis, we reconsider the partial differential 
equations ( |3.15| ) this time for small values of the initial couplings. Expanding 
F and G as 

F(<p;t) = fi 2 (t) + \(t)<p 3 + ... 

G{ V -t) = g(t)(p + ... (3.24) 

and retaining only the lowest order terms, the differential equations reduce 
to 

l 6vr 2^) = _i8A 2 (t)+% 4 (t) 

16tt 2 ^ = -2g\t) (3.25) 
at 

which are precisely the one-loop perturbative renormalization group equa- 
tions in the approximation in which the anomalous dimensions are neglected. 
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In obtaining these equations, we have assumed that the parameter /x 2 (0) has 
been tuned to be close to its critical value, /j% rit = — A(0)(t> 2 /A 2 ), which is 
proportional to A(0) and hence small. Here we denoted the <p vacuum ex- 
pectation value by v/A. The analytic solution to these equations is readily 
secured as 



1 



9 2 (t) 
X(t) 



+ 



4vr 2 



t 



a(t) 



where 



a(t) 



37(a(t) - 1) - (a(t) + 1) 
37+1+40 

4 A(0) 



37 



(3.26) 



(3.27) 



Extrapolating these solutions well outside their domain of validity, we display 
these perturbative results for the running couplings as the dashed lines in 
Figs. 4-5 using the same initial inputs as used in the numerical solution 
to the nonperturbative locally approximated Wilson renormalization group 
equation. 

As can be seen in the g(0) = 5, A(0) = 1 case, the perturbative renormal- 
ization group running of X(t) and g(t) exhibits the same qualitative features 
as the nonperturbative evolution, but are quantitatively quite different. The 
large Yukawa coupling again initially causes A(i) to increase. This time, the 
growth is quite dramatic, as depicted by the dashed line curve in Fig. 4, 
since the higher dimensional operators or even higher loop self coupling ef- 
fects which provide the additional damping are now absent. This rise in X(t) 
continues until g(t) has decreased and A(£) has increased sufficiently so that 
g 2 (t) = |A(t). At this point the slope of the curve is zero. Subsequently, the 



pure scalar self coupling term in Eq. (|3.25|) drives A(t) slowly towards zero. 
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Figure 4: Comparison of the numerically integrated locally approximated 
Wilson (solid line) and 1-loop perturbative (dashed line) renormalization 
group equation solutions for the running of the scalar self coupling for various 
initial Yukawa couplings and cutoff to scalar vacuum expectation value ratios. 

For the case of a smaller initial Yukawa coupling constant and larger 
initial scalar self coupling, g(0) = 1, A(0) = 5, g(t) is seen in Fig. 5 to 
be completely described by its perturbative renormalization group evolution, 
while Fig. 4 shows X(t) monotonically decreasing from its initial value since 
g 2 (0) < y 1+^2(0) ■ This is again in crude qualitative agreement with the exact 
solution, although the decrease in X(t) in the exact solution is far more rapid 
due to the effects of the induced higher dimensional operators (cf. Fig. 1). 
On the other hand, for the larger cutoff, A = 38. 6i>, the effects of the higher 
dimensional operators are far less significant and the t-evolution of X(t) runs 
essentially as in 1-loop renormalization group improved perturbation theory. 
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Figure 5: Comparison of the numerically integrated locally approximated 
Wilson (solid line) and 1-loop perturbative (dashed line) renormalization 
group equation solutions for the running of the Yukawa coupling for various 
initial scalar self couplings and cutoff to scalar vacuum expectation value 
ratios. 
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4 Scalar and Fermion Mass Bounds 



The numerical solution to the locally approximated Wilson renormalization 
group equation can also be used to obtain bounds on the domain of allowed 
scalar and fermion masses. Working in the broken 75-symmetry phase, we 
define the f and ip masses as the zero momentum value of the one particle 
irreducible scalar and fermion two point Schwinger functions, r( 2 ' 0,0 )(p) and 
r^ ' 1 ' 1 )^), respectively, so that 

^ r (2,o,o )(p = 0) (41) 

^- = \trT^(p = 0). (4.2) 

Since Z = 1 in the local action approximation we are employing, these masses 
are also the location of the poles of the propagators when analytically con- 
tinued to Minkowski space. 

The short distance contributions to the Schwinger functions are evalu- 
ated by numerically integrating the Wilson renormalization group equation 
from the initial cutoff A, where the generalized potential is specified, down 
to a lower momentum scale, A* = e~**A, below which the dynamics can be 
accurately described perturbatively. The contributions to the zero momen- 
tum Schwinger functions from degrees of freedom with momentum below A* 
(infrared contributions) are then included as (one-loop) perturbative correc- 
tions to the generalized potential U(ip, a; t*) obtained from the numerical 
integration. In this sense, the Wilson renormalization group equation pro- 
vides a systematic procedure of generating the effective Lagrangian valid on 
one distance scale from that of a more finely grained distance scale. Integrat- 
ing all the way down to zero momentum (t — > 00), the effective generalized 
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potential is secured as (cf. Eqs. ( |3.21| - |3T23| )) as 



U eff (ip, 



where 

V eff (<p) 



V eff (<p) + aG eff (<p), 



str In 



d 2 U(^,(r;t*) 



(4.3) 



16tt 2 



64tt 2 



(l-^(^;r))ln + 



+^(^5 f ) + Vlfa; t*)ln V^fa; t* ) 



(1 - G%; f )) In (1 + G 2 (<p; 1*)) + G 2 (^; f ) 



+G 4 (^;t*) In G 2 (^t*; 



(4.4) 



32tt 5 



G w (<^**) + 



2G(^t*)Gl( V -t*) 



l+VUw?) In 



(\/ w (^t*)-G 2 (^t*))J 



G((f, t*)Gl((p; t* 



16tt 2 {y<p<p{<p',t*) — G 2 ((p;t*)) 



1 + G 2 {<p;t*) In 



(4.5) 



To proceed with the numerical integration of Eq. Q3.15 ), we first determine 
the critical value of the mass parameter, fi 2 rit , for all possible choices of the 
initial coupling constants A(0) and g{0). For each particular choice of initial 
A(0) and g(0) parameters, n 2 crit is defined so as to separate the discrete 75- 
symmetric phase (/i 2 (0) > /i 2 rit ) from the spontaneously broken 75- symmetry 
phase (/i 2 (0) < fi 2 rit ). The broken symmetry ground state is determined by 
the location of the zeros of F((p;t). The value of [i 2 crit is the maximum of 
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/i 2 (0), for a given A(0) and g(0), which results in a nontrivial zero of F(ip; t) 
as t increases into the infrared. This is determined by evaluating F(<p; t) for 
/i 2 (0) well into the broken phase and then increasing /i 2 (0) until the vacuum 
value decreases, rather than increases, as the model is evolved in t. The 
transition value of yu 2 (0) delineating the increasing and decreasing evolution 
of the vacuum value is /i 2 rii . Performing this analysis, one also finds that the 
model has only a trivial (Gaussian) infrared fixed point in accord with the 
conclusions based on lattice simulations. ^ 

Once the critical line is determined, a choice of initial parameters A(0), 
g(0) and /i 2 (0) < /i 2 rit is made and the t dependence of the generalized poten- 
tial terms F(<p; t) and G(<p; t) is extracted. This nonperturbative evolution is 
continued until the the scale t* is reached at which point the dynamics of the 
system can be accurately described by perturbation theory and the general- 
ized potential is given via Eqs. ( |3.21[ - |3l23| ). The scalar and fermion masses 
are then extracted using the (one-loop) perturbative expressions for the cor- 
responding zero momentum Schwinger functions. Thus the scalar mass is 
given by the curvature of the effective potential at its minimum via 

= V#>(~) (4.6) 

and the fermion mass is simply the effective Yukawa coupling constant also 
evaluated at the minimum of the effective potential: 

= G^(j). (4.7) 
The vacuum expectation value, v/A, is given, through one loop, by the 
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location of the minimum of the effective potential as 



v; e// (^) = o. 



(4.8) 



Letting Vq denote the vacuum expectation value of the the tree level potential 
V((p; t*), then the one loop shift in expectation value, 5v , is just 



5v 
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4tt 2 VI, 
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;i + g 2 ) 



(4.9) 



tp=v /A, t=t* 

Using these results in Eqs. ( f4.4| - f4.7| ), the scalar and fermion masses are 
secured in terms of the effective potential and effective generalized Yukawa 
coupling evaluated at t* and the tree vacuum value vq as: 
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The value of t* is ascertained by the desired numerical accuracy of these 
masses determinations. As the system evolves into the perturbative regime, 
the mass values obtained remain stable as t* is varied. We choose for t* that 
particular value such that the computed mass is stable to within an error the 
size of the t-grid spacing: 10~ 5 < 5t < 1CT 4 . 

Using the above procedure, the ratios M s /A, mf/A and v/A are deter- 
mined as the initial couplings A(0) and g(0) are varied. Taking appropriate 
products of these ratios, we plot in Fig. 6 the allowed M s /v values as a 
function of A/M s for various different choices of A(0) and g(0). Note that 
the plotted data is further restricted to satisfy the physically imposed con- 
straint that A/M s > 1. This follows since no masses are allowed to become 
larger than the initial UV cutoff A. An upper bound on the M s /v values as 
a function of A/v is secured by noting that a dynamical envelope is formed 
as the initial A(0) and g(0) couplings are varied. That is, for fixed A/v, the 
different M s /v values obtained appear to converge to an upper limit, at least 
for initial Yukawa couplings < 10. i 

Focusing on the dependence of M s /v on the initial Yukawa coupling, we 
note that for A(0) = 5, a change in g(0) from to 7 leads to roughly a 
15-20% increase in M s /v. On the other hand, for the larger initial value of 
A(0) = 15, the same change in g(0) from to 7 results in only an 6-10% 
increase in M s /v. Thus the effect of the initial Yukawa coupling on M s /v 

3 Due to limitations in the presently used numerical integration technique, we do not 
establish the dynamical envelope for larger values of the initial Yukawa coupling g(0). 
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Figure 6: M s /v as a function of A/M s for different initial A(0) and g(0) 
values. 
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diminishes with increasing initial A(0). 

For fixed A(0) and g(0), note that all the curves in Fig. 6 except the 
one denoted with the crosses (x), A/v increases as A/M s increases. The 
contrary situation persists for the curve with A(0) = 0, g(0) = 7 where the 
maximum scalar mass also corresponds to the largest A/v. Furthermore, 
the curve corresponding to A(0) = 5 and g(0) = 7 has the behavior that 
M s /v initially increases with A/v before reaching its maximum and then 
subsequently deceases. Except for these two curves, which are characterized 
by having g(0) > A(0), all the other curves exhibit a maximum scalar mass 
when A/v ratio is at its minimum. The reason for the behavior of these two 
curves is that the initially large g(0) value requires a longer evolution time 
(or equivalently a larger initial cutoff) to indirectly feed into an increase in 
the scalar mass. 

In addition, we can extract the allowed domain of scalar and fermion 
masses for a fixed value of the ultraviolet cutoff to scalar vacuum expec- 
tation value ratio A/v. Figure 7 is just such a plot ' 13 ' corresponding to 
the specific choice A/v = 5. The allowed values are those lying interior to 
the various boundaries displayed. We now discuss the origin of each of the 
boundary restrictions in turn. The upper portion of the boundary, denoted 
by the cross (x) marks, is the scalar mass triviality bound discussed above. 
Note that the absolute scalar mass upper bound, for this value of A/v , occurs 
for the larger fermion masses and is simply given by M s = A = 5v. Next, 
the square markings (□) correspond to the triviality bounds on rrif. These 
points are obtained by finding the largest g(0) values for a given A(0) which 
produce an envelope in the rrif/v versus A/rrif plane. Finally, by demand- 
ing that the scalar potential, V(ip;t), be bounded from below on all scales, 
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Figure 7: Allowed range of M s /v and rrif/v values for A/v — 5. 



we further restrict the allowed masses. This vacuum stability requirement, 
which corresponds to the condition that lim^^^F^ip; t) > for all t, is 
implemented by setting A(0) = and allowing g(0) to vary. This results 
in the points on the boundary in Fig. 7 represented by the diamonds (O). 
Note that for larger fermion masses, the fermion mass triviality condition 
provides a more stringent bound than that of vacuum stability. This is not 
the case for very large A/v ratios where the vacuum stability constraint is 
always more stringent 1 13 1 than that due to the fermion triviality. 
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5 Conclusions 



The Wilson renormalization group equation provides a systematic nonper- 
turbative tool for constructing the action at a particular distance scale given 
the form of the action at a more finely grained scale. The resulant action con- 
tains operators of all mass dimension. In particular, if one defines the form 
of the action at some UV cutoff scale A, then for slightly lower scales, the ir- 
relevant operators (those of mass dimension greater than the dimensionality 
of space-time) can play a highly nontrivial role in the ensuing dynamics. 

We constructed the Wilson renormalization group equation for a theory 
involving a scalar and a fermion and numerically solved the equation in a 
local approximation. The resulting renormalization group trajectories were 
analyzed and we established triviality and vacuum stability mass bounds on 
both the scalar and fermion degrees of freedom. Clearly, it is desirable to 
extend this study beyond the local approximation to include the anomalous 
dimensions (nontrivial wavefunction renormalizations) for the scalar and fer- 
mion as well as higher dimension derivative couplings. At the present time, 
the complications limiting such an extension are purely of a numerical na- 
ture and we plan to address these issues in future studies. In addition, the 
analysis can also be enlarged to allow for the more realistic case of scalars 
and chiral fermions transforming nontrivially under an internal non-Abelian 
Lie group and also via the inclusion of gauge bosons. This later step can 
prove technically challenging since the hard cutoff we employ here breaks 
the manifest gauge invariance. Finally, the model can also be made super- 
symmetric and an analysis similar to that performed here can yield SUSY 
mass triviality bounds. 
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Appendix: 

Supermatrices, Supertraces and Superdeterminants 

When discussing theories containing both bosons and fermions, it often 
proves convenient to combine terms possessing different Grassmann algebra 
character into a single supermatrix structure. In this appendix, we review 
some of the properties of the supermatrix 

/ M BB M BF 

M = 

V M FB Mpp 

where the square submatrices M BB and Mff contain only even elements 
of a Grassmann algebra while M BF and M FB contain odd elements of the 
Grassmann algebra. In particular, M BB and M F f are so defined such that the 
components of M BB are composed of some ordinary c- numbers (Grassmann 
number zero) as well as possibly some higher even elements of the algebra, 
while the components of Mff are of Grassmann number two or higher. 
We define the supertrace (str) of the supermatrix M via 

str M = tr M BB - tr M FF , (A.2) 

where tr denotes the ordinary trace. This definition is chosen so that the 
familiar trace property 

str (MN) = str (JVM) (A.3) 

is still guaranteed to hold. Here the multiplication law for supermatrices is 
defined precisely the same as for ordinary matrices. Using the definition of 
the str operation, we next define the superdeterminant (sdet) by 

sdet M = e str ln M . (A.4) 




36 



Using the str property (A. 3), it is straightforward to establish that 



<5(sdet M) = sdet M ■ str (M _1 5M) 

and that 

sdet (MAO = sdet M ■ sdet N ■ 
Now consider decomposing the supermatrix M as 

(M BB 0\ (\ M BB M BF \ 



M 



M. 



FB 







N. 



FF 



where 



N FF = M FF - M F B Mg B M BF . 
From property (A. 6), it follows that 

1 Mz l B M BF - 
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(A.5) 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



FF 



Because each of the supermatrices appearing on the right hand side in the 
decomposition ( |A.7| ) has zeroes in either the upper right or lower left off 



diagonal blocks, their respective sdet can be directly evaluated using the 
definition ( |A.4| ) and the str property ( |A.2| ). For example, 



sdet 



M BB 0' 
. M FB 1 , 



exp < str In 



1 
1 

exp|strE^ =1 ^ 
exp {tr En=i ^ 



M BB -1 0' 
M FB 0. 



M BB -1 0' 

V M FB 0, 
(M BB - 1)"} 



exp {tr \iaM BB } 
det M BB , 



(A.10) 
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where det is the ordinary determinant. Similarly, one finds 

/ 1 MzlM BF \ 
sdet = exp{(-)tr lnN FF } 



N, 



FF 



(det Nff)- 1 . (A.ll) 



Here the negative sign appearing in the exponential, which in turn leads to 
the inverse determinant, is a direct consequence of the supertrace definition, 
Eq. (|A.2|) . Combining terms, we secure 

det M BB / » -„\ 

sdet M = - -fl . A. 12 

det N F p 

Alternately, we could have chosen to decompose M as 

(N BB M bf Mff\ ( 1 \ 
M=[ J • I , A. 13 

V 1 J \M F B M ff J 

with 

-Abb = M BB - M BF Mp F M FB . (A. 14) 

So doing, one finds 

det M = P Nbb . (A.15) 
det M FF 
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